Hadronization and Hadronic Freeze-Out in Relativistic Nuclear Collisions 
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We analyze hadrochemical freeze-out in central Pb+Pb collisions at CERN SPS energies, employ- 
ing the hybrid version of the Ultrarelativistic Quantum Molecular Dynamics model, which describes 
the transition from a hydrodynamic stage to hadrons by the Cooper-Frye mechanism, and matches 
to a final hadron-resonance cascade. We fit the results both before and after the cascade stage 
using the statistical model, to assess the effect of the cascade phase. We observe a strong effect on 
antibaryon yields except anti-fi, resulting in a shift in T and \ib of the freeze-out curve. We discuss 
indications of a similar effect in SPS and RHIC data, and propose a method to recover the bulk 
hadron freeze-out conditions. 
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I. INTRODUCTION 

Hadron production in relativistic A+A collisions is 
supposed, since Bevalac times [l|, to proceed via two 
separate stages. The first, "hadrochemical" freeze-out 
fixes the bulk hadronic yields per species which are con- 
served throughout the subsequent hadron-resonance cas- 
cade expansion. At its end, "kinetic freeze-out" delivers 
the eventually observed bulk properties such as pt spec- 
tra, HBT correlations, collective flow, etc.. Most remark- 
ably, the hadronic yield distributions over species are un- 
derstood to resemble grand canonical statistical Gibbs 
equilibrium ensembles |2|, H[ , from AGS up to RHIC and 
LHC energies. The two most relevant parameters, tem- 
perature T and baryochemical potential (1b thus capture 
a snapshot of the system dynamical evolution, taken at 
the instant of hadrochemical freeze-out. 

In relativistic A+A collisions the thus determined T 
increases monotonically with -Js NN , saturating at about 
170 MeV while /is approaches zero. Systematic statisti- 
cal model (SM) analysis reveals the "freeze-out curve" Q 
in the T, \ib plane, in which we usually also represent the 
conjectured plot of the phase diagram of QCD matter. 

Such a plot is given in Fig. Q] It shows two principal 
lines, firstly a parton-hadron coexistence boundary line, 
inferred from lattice QCD [j| at low /is , and from chiral 
restoration theory Q at high [Ib- And, second, the SM 
freeze-out curve. Remarkably, the lines merge toward 
T = 170 MeV, /is = 0. The freeze-out curve thus lo- 
cates the QCD hadronization transition temperature T c : 
hadronization seems to coincide with hadronic freeze-out, 
here [7j. Equally remarkable, however, the two lines dis- 
entangle with increasing /is, becoming spaced by about 
30 MeV temperature difference toward [Ib = 500 MeV 



> 



200 



quark gluon plasma 



100 



E 



chemical frcczc-out 
A SIS, AGS 
■ SPS 
RHIC 



hadron gas 



quark- 
yonic 
matter 



M 



500 



1000 
MMeV) 



"now at Yale University, New Haven, CT, USA. 



FIG. 1: (Color online) Sketch of the QCD phase diagram, 
including the hadronic freeze-out curve (see text). 



which corresponds to ^Js NN — 5 GeV in A+A collisions. 

What are we freezing out from, here? If hadronization 
coincides with the parton-hadron transition line at small 
/j,b, and if hadronization involves a statistical equilibrium 
of hadron species Q , one could imagine that the suc- 
cess of the statistical model in reproducing hadronic mul- 
tiplicities in the higher fj,s region indicates that the chem- 
ical freeze-out line marks the transition to the hadron 
world from an alternate phase. Indeed, the possible exis- 
tence of an alternative QCD phase, so-called quarkyonic 
matter llOf has been put forward based on similar ar- 
guments [11] . Indicated in Fig. Q] is a scenario in which 
the freeze-out curve is identified, tentatively, with a hy- 
pothetical quarkyonic matter phase boundary. 

Before embarking on this idea a different possible sit- 
uation needs to be addressed. Taking for granted that 
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the hadron-resonance phase is indeed created directly at 
the coexistence curve it might be conceivable that an 
expansive hadron/resonance evolution stage, setting in 
at hadronization temperature and baryochemical poten- 
tial T c and hb,c, cools down the population maintain- 
ing chemical equilibrium until chemical freeze-out occurs 
by mere dilution (the inelastic mean free path becoming 
larger than the system size), at lower T and higher \xb-, 
thus defining the freeze-out curve. 

In this paper we test the latter scenario. We employ 
the framework of the Ultrarelativistic Quantum Molec- 
ular Dynamics (UrQMD) microscopic transport model. 
Its recent hybrid version [12J features a 3+1-dimensional 
hydrodynamic expansion during the high density stage, 
terminated by the Cooper-Frye mechanism once the en- 
ergy density of flow cells falls below a "critical" energy 
density, assumed to be e c = 0.8 GeV/fm 3 . This crite- 
rion resembles the parton-hadron transition line of Fig. [T] 
The hadron/resonance population can be examined, ei- 
ther, by terminating the model evolution at this stage, 
emitting into vacuum, and fitting the yield distribution 
by the grand canonical statistical model [.13]. Or, alter- 
natively, the UrQMD hadron/resonance cascade expan- 
sion stage is attached, as an "afterburner" . The outcome 
is again fitted by the SM. A preliminary report on this 
work was published in Ref. We shall show that 

the afterburner does not cool the system in equilibrium 
(to start from T c and arrive at a new equilibrium Tf on 
the freeze-out curve). However, it is shown to selectively 
absorb antiprotons and antihyperons except anti-Os at 
SPS energies, while leaving the other yields essentially 
unchanged. A distorted yield distribution results, as it 
was already shown by Bass and Dumitru 15j . We shall 



explore a method to circumvent such final state effects, 
recovering the true hadronic freeze-out curve. 



II. URQMD 
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FIG. 2: (Color online) UrQMD calculations for hadron mul- 
tiplicities in central Pb+Pb collisions at (a) 40A and (b) 
158A GcV, plotting the result before the hadronic cascade vs. 
the result after. Antibaryons are highlighted by red markers. 



Figure [5] shows the effect of the final UrQMD cascade 
stage, in a plot of hadron multiplicities directly after the 
hydro stage, vs. the multiplicities at the end of the cas- 
cade. We illustrate these conditions for the 5% most cen- 
tral Pb+Pb collision multiplicities, at the SPS energies 
40 and 158A GeV (^s^ = 8.6 and 17.3 GeV). We see 
the bulk hadrons unaffected by the afterburner, includ- 
ing the S and anti-f2. Whereas the other antibaryons, 
p, A and 2, are significantly, and selectively suppressed. 
We quantify this suppression in Fig. 3 showing the var- 
ious hadron multiplicity ratios, of full UrQMD calcula- 
tions ("Hydro+ Afterburner"), relative to the output af- 
ter Cooper-Frye transition ("Hydro"), for the two cases 
illustrated in Fig. [2] The p, A and 5 multiplicities show 
reductions ranging from 45 to 25%, whereas the other 
hadron yields stay essentially constant. Similar results 
are reported in ref. 15[ where it was shown, further- 



more, that at top RHIC and LHC energies the selective 
absorption of these antibaryon species turns into a gen- 



eral suppression of both the corresponding baryon and 
antibaryon species. This should be a consequence of the 
approximate particle-antiparticle symmetry occurring at 
such high energies, where \ib approaches zero. 

The mechanism of this suppression of antibaryons 
(at SPS energies), in the course of the UrQMD 
hadron/resonance cascade phase, appears to be the ex- 
citation to higher resonance states which increases the 
inelastic (annihilation) interaction rate, over and above 
"direct" antibaryon-baryon annihilation. The resonance 
mechanism should not be active for the f2 hyperons which 
lack such low lying excitations. 

Figure [3] also illustrates similar UrQMD results ob- 
tained for semi-central Pb+Pb collisions in the 35 to 45 
percent centrality bin where (N part ) is about 90. Qualita- 
tively similar suppression occurs, albeit of reduced mag- 
nitude owing to the shorter duration of the final cascade 
stage. We shall return in Sec. HVI to experimental conse- 
quences of this annihilation increase toward central col- 
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FIG. 3: The ratios of hadron multiplicities after the UrQMD 
cascade stage, to the multiplicities directly after hadroniza- 
tion (after the hydro-stage). We illustrate central Pb+Pb 
collisions at (a) 40 and (c) 158A GeV, as well as the corre- 
sponding semi-central cases for {N pal t) = 90 (b), (d). 



lisions. 



III. STATISTICAL MODEL ANALYSIS AND 
THE FREEZE-OUT CURVE 

Turning to statistical model analysis of the UrQMD 
results we consider the total hadron multiplicities in cen- 
tral Pb+Pb collisions at the SPS energies 20, 30, 40, 80 
and 158A GeV (y^7 = 6.3, 7.6, 8.7, 12.3 and 17.3 GeV). 
They are chosen to obtain a sufficient coverage of the 
freeze-out curve illustrated in Fig. [1] and to resemble the 
positions of the data points gathered by NA49 [THH?]]- In 
order to obtain results under statistical conditions similar 
to those prevailing in the analysis of the corresponding 
data we attach to the hadron multiplicity results from 
UrQMD the corresponding total error resulting from the 
addition in quadrature of the statistical and systematic 
errors reported by N A49 [l|| ■ We employ the grand 



canonical version of the statistical model as described in 
refs. @, [Hj]. In this version, the statistical model is sup- 
plemented with a free parameter 75 suppressing the pro- 
duction of hadrons containing n s valence strange quarks 
according to 7^ . In the SPS energy region, as well as in 
the peripheral collisions at RHIC energy fTsL HH this pa- 
rameter turns out to be less than one [2|, Il3| . It has been 
observed that this parameter can be explained, at least 
in the RHIC energy region, with the superposition of two 
sources (the core-corona model): single nucleon-nucleon 
(NN) collisions and a fully equilibrated source (2(j. For 
all hydro+UrQMD calculations, this parameter is set to 
1 and no core-corona effect is included. 

Figure |4] illustrates the fits to the hybrid UrQMD re- 
sults by the statistical model, choosing the 158 A GeV 
case as an example. The afterburner stage indeed shifts 
(T,fi B ) considerably, from (160 MeV, 246 MeV) to 
(151 MeV, 277 MeV). However, note the dramatic de- 
crease of fit quality, from 4.7 to 28.2 in chisquare. Sim- 
ilarly, the results for the energy 40A GeV are shown in 
Fig 03 exhibiting the same shift conditions. The effect 
of the afterburner is, thus, not an in-equilibrium cooling 
but rather a distortion of the hadron yield distribution 
in the antibaryon sector, away from equilibrium - as we 
could guess from Figs. [2] and El already. The parameter 
7s is very close to 1 throughout. 

The idea arises to exclude p, A and 2 from the SM 
fit. Figsures [6] and [7] show examples, again at 158 and 
40A GeV. No cooling occurs. The fit to the afterburner 
output (which features a reasonable chisquare) now ig- 
nores the far off-diagonal p, A and 5 entries. A sum- 
mary of the results, obtained for all energies, is given in 
Fig. [8] The (T,hb) positions of all analyzed cases are 
represented in a phase diagram that also exhibits the av- 
erage freeze-out curve resulting from previous statistical 
model analysis of experimental data |4j. We see that 
the series of entries resulting from the full UrQMD cal- 
culation, including the final hadron/resonance cascade 
stage, follows rather closely the freeze-out curve. The 
points obtained without the cascade stage are generally 
well above, and also shifted to lower fig. Not surpris- 
ingly, they approach the lattice parton-hadron boundary 
shown in Fig. [TJ Most significantly, however, the points 
from full UrQMD but with omission of the p, A and S 
multiplicities are also following the latter behavior: they 
universally exhibit no significant "cooling" , with the ex- 
ception, perhaps, of the high /is region. 

We conclude that the hadron/resonance cascade as 
modeled in the microscopic dynamics of UrQMD can not 
transport an initially established hadrochemical equilib- 
rium from the phase coexistence line of Fig. [T] downward 
to the freeze-out line. However, it distorts the hadron 
yield distribution which leads to a downward shift of 
the freeze-out parameters derived from statistical model 
analysis, albeit at the cost of a rather unsatisfactory fit 
quality. Far better fits are obtained omitting the three 
antibaryon species, which shifts the freeze-out curve up- 
ward. We note, however, that it will be of key interest 
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FIG. 4: (Color online) Statistical model fit to UrQMD re- 
sults for hadron multiplicities in Pb+Pb central collisions at 
158A GeV. (a) employs UrQMD terminated directly after the 
hydro stage whereas (b) refers to multiplicities obtained after 
the final hadron/resonance cascade phase of UrQMD. 



FIG. 5: (Color online) Same as Fig.[4]but for Pb+Pb collisions 
at 40 A GeV. 



to extend the present study toward smaller /is. Will 
the three curves in Fig. [8] converge, more closely, toward 
[iB = 0? This is a subject for ongoing study. How- 
ever, both former (l5j and concurrent 21 1 investigations 
with a hybrid UrQMD model conclude on significant af- 
terburner effects, even up to LHC energies. Along with 
forthcoming data (see next section) the method proposed 
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FIG. 6: (Color online) Same as in Fig. [4] but excluding p, A 
and 3 from the statistical model fit. Note that the residual 
plots do nevertheless also illustrate the entries resulting for 
these antibaryons. 
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FIG. 7: (Color online) Same as Fig. ©for Pb+Pb at 40A GeV. 



IV. A LOOK AT THE DATA 



We note that the properties reported above may be 
a consequence of the specific hadron transport model 
employed in the UrQMD approach, which implements 
here may thus lead to a certain, general upward shift of a cascade of binary hadron/resonance collisions, com- 
the freeze-out curve at low /is- bined with string excitation in such collisions which leads, 
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FIG. 8: (Color online) Summary of the results obtained for 
central Pb+Pb at five energies in the SPS range (^/s NN — 6.3, 
7.6, 8.7, 12.3 and 17.3 GeV). T and /is are shown from statis- 
tical model fits, at each energy, to the three UrQMD config- 
urations illustrated in Figs. [4] to [71 The dashed line indicates 
the "empirical freeze-out curve" [4J. The fit to hadron mul- 
tiplicities right after Cooper-Frye transition from the hydro 
stage (squares) results in higher temperatures. The after- 
burner stage of hadron/resonance rescattering cools the sys- 
tem (black triangles). The original freeze-out conditions are 
essentially recovered by the modified fit (red triangles) as de- 
scribed in the text. 



TABLE I: Main fit to hadron multiplicities (in full phase 
space) measured by NA49 in central Pb+Pb collisions at 
\fs NN = 17.3 GeV 0,113] (Fig EI (a)). Errors within brack- 
ets are the (realistic) parameter errors, obtained by rescaling 
the original fit errors by \/ y 2 /DOF according to the method 
suggested by the PDG [2fJ. 

Fit Parameters 

T (MeV) 150.5±1.7 (3.1) 
7s 0.915±0.034 (0.063) 

fin (MeV) 266±6 (11) 
X 2 /DOF 26.9/10 



Hadron Measured 


Fitted 


B 


362±8 


367 


7T+ 


619±35 


528 


7T~ 


639±35 


560 


K+ 


103±7 


105 


K" 


51.9±3.6 


60.5 


k" 


81±4 


82 


A 


48.5±8.6 


51.1 


t 


8.46±0.50 


8.27 




4.40±0.64 


4.80 


P 


4.23±0.35 


4.69 


A 


3.32±0.34 


3.17 


e+ 


0.710±0.098 


0.592 


n 


0.59±0.11 


0.49 


n 


0.26±0.07 


0.13 



in part, to multiparticle decays. However, the detailed- 
balance counterparts of the latter processes are not in- 
cluded. It has been argued repeatedly d, that such 
multi-hadron collisions could significantly decrease the 
effective chemical relaxation constants governing antihy- 
peron densities and their approach to chemical equilib- 
rium. However, in these approaches only the direct vicin- 
ity of the critical temperature, T c — 160 — 170 MeV, has 
been addressed. Thus this line of argument does, in fact, 
refer to the parton-hadron phase transformation itself. It 
attempts to recast the hypothesis of equilibrium estab- 
lishment via the quantum mechanical mechanism of the 
phase transition 0, d, HH m terms of a quasi-classical mi- 
croscopic dynamics of a mixed phase, consisting of Hage- 
dorn resonances or string excitation modes and hadrons 
plus hadronic resonances, in some sense a quasi-classical 
microscopic mixed phase, and phase transition model. 
Such models, however, work in the immediate vicinity of 
T c only, i.e. in application to RHIC and LHC energies. 
In the present UrQMD study we address, specifically, the 
possibility of a hadrochemical equilibrium freeze-out at 
significantly lower temperature, and high baryochemical 
potential, where the freeze-out curve disentangles from 
the parton-hadron phase boundary. It is our observa- 
tion that the UrQMD model dynamics can not trans- 
port the system, from initially imprinted equilibrium by 
the Cooper-Frye mechanism, while maintaining chemical 
equilibrium adjusted dynamically to the rapidly falling 



temperature during expansion. This is in line with a 
former study [24j of UrQMD equilibrium features in the 
same energy domain. Putting the system into a box, with 
constant energy density, it was shown that chemical equi- 
libration requires upward of 25 fm/c, even under constant 
energy density conditions, not available during cascade 
expansion. From that study we would conclude that, 
conversely, an equilibrium initially established would not 
be substantially distorted during the afterburner phase. 

In the (T,(1b) domain treated here, the transport 
model employed in UrQMD may appear to be realistic. 
The idea thus arises to look for effects, similar to the 
findings summarized in Fig. [51 in real data. We recall 
from previous statistical model analysis of the SPS to- 
tal multiplicity data @, [l3[ that the obtained chisquare 
values tended to be rather high. We repeat here the SM 
analysis of the NA49 up-to-date data set (see Table [T| for 
central Pb+Pb collisions at top SPS energy, 158 A GeV 
= 17.3 GeV). With the updated data sample by 
NA49, the above trend is confirmed: the fit quality is 
worse than in previous analyses as seen in Fig. [S] (a) and 
Table [TJ This might be an indication that the inelastic 
rescattering stage must be taken into account if one aims 
at matching the accuracy of the experimental measure- 
ments. We have not repeated the analysis at lower SPS 
energies because the coverage of the hyperon sector was 
incomplete in NA49 but we note that new high precision 
results at y/s NN = 7.7 GeV are forthcoming from STAR, 
gathered in the recent low energy runs at RHIC [26[. 
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FIG. 9: (Color online) Statistical model fit to the NA49 
data [IB, 03] for total hadron multiplicities in central Pb+Pb 
collisions at 158^4 GeV. (a) full data set included in the fit. 
(b) fit with excluded p, A and 3 entries (note that, neverthe- 
less, the resulting residuals are illustrated). 



We show in Fig. [9] a comparison of a SM fit to the full 
set of 158 A GeV NA49 data [H, [T3| (Fig. E] (a)), with a 
fit where p, A and S were excluded (Fig. [9] (b)). Indeed, 
in accordance with the corresponding UrQMD exercise 
shown in Fig. [4l the resulting SM parameters move up 
in the temperature T by about 25 MeV, in the latter 
case, along with an improved chisquare. However, this 
large shift is accompanied by a significant decrease of 75, 
which makes the conclusion of a major role of afterburn- 
ing in antibaryon absorption questionable, because such 
a shift is not observed in the simulated data, as has been 
mentioned. In order to put our conclusion on a firmer 
ground, we have repeated the fit by fixing 75 at a pre-set 
value of 0.85 (see e.g. Refs. (IE!)' 

in both cases. The re- 
sults are shown in Fig. [10] (top and bottom panels). The 
shift in temperature when excluding antibaryons (except 
f2) is smaller than in the fits shown in Fig. [9j but still 
significant. 

We show in Fig. QT] that the differences between the 
data and the statistical model fit, shown in Fig. [TU] (b), 
exhibit a striking resemblance to the corresponding 
"hadron survival" fractions, obtained from the UrQMD 
study illustrated in Fig. [3] Selective suppression of p, A 
and 5 thus appears to be a property of the data. 

A further indication of selective antibaryon absorp- 
tion can be found in STAR Au+Au minimum bias data 
at 7.7 Ge V J26| . and in NA49 minimum bias data at 
17.3 GeV jl7l.l27|. The kaon, A and S multiplicities per 
participant nucleon increase smoothly toward central col- 
lisions, a property that can be well explained by the de- 
creasing relative corona contribution 20, 28, 29]. On the 
contrary, the p, A and H yields per participant stay con- 
stant with centrality in these collisions. One might thus 
infer on a further centrality dependent effect that works 
selectively in the latter cases. We recall the observation 
made in Fig. [3] that the antibaryon absorption increases 
with centrality. This could counterbalance the yield in- 
crease per participant nucleon, exhibited by the baryons. 

At the much higher top RHIC, and LHC energies, the 
selective suppression of the antibaryon yields (except for 
the hyperon) in the hadron/resonance cascade stage, as 
implied in our results at SPS energies, should give way to 
a suppression symmetric in baryons and antibaryons (as 
was already shown in Ref. [Hj]), due to the approximate 
particle-antiparticlc symmetry prevailing here. Measure- 
ments of proton and antiproton midrapidity yields at 
LHC seem to confirm this [3(i| . It is a likely possibility 
that the exclusion of those baryons from statistical model 
fits at higher energy will also lead to an improvement of 
the fit quality and to a slight increase of the estimated 
chemical freeze-out temperature values. 



CONCLUSIONS 



We have systematically investigated the effect of the 
hadron-resonance expansion phase that follows hadro- 
chemical freeze-out, employing the hybrid version of the 
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FIG. 10: (Color online) Statistical model fit to the NA49 
data p3, E3] for total hadron multiplicities in central Pb+Pb 
collisions at 158A GeVwith 7s = 0.85 fixed, (a) full data 
set included in the fit. (b) fit with excluded p, A and 3 
entries (note that, nevertheless, the resulting residuals are 
illustrated) . 
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FIG. 11: Hadron multiplicity ratios data/SM implied by the 
statistical model fit to the data set excluding the three an- 
tibaryons, shown in Fig. (b). 



UrQMD transport model, which contains a transition 
process (from an intermediary hydrodynamic evolution), 
as described by the Cooper-Frye mechanism. We con- 
sider the hadron multiplicities obtained with, and with- 
out the cascade phase, analyzing them with the statistical 
model. We find that the cascade expansion as modeled in 
the microscopic dynamics of UrQMD does indeed modify 
the hadronic yield distribution, thus significantly affect- 
ing the position of the freeze-out curve obtained from 
the SM. A downward shift occurs in the (T,/j,b) plane 
that, at first sight, appears to bridge the gap between 
the position of the lattice hadronization curve, and the 
well known statistical model freeze-out curve, at least at 
the investigated SPS energies. However, a closer anal- 
ysis of the SM results shows that the observed shift is 
universally caused by a selective suppression of the an- 
tibaryons (except f2) during the cascade expansion stage. 
The other hadron yields stay very closely constant. What 
we observe is, thus, not a cooling of the hadron distri- 
bution but a distortion, and in fact it is accompanied 
by a systematic decrease of SM fit quality. If one omits 
these antibaryons from the SM fits most of the downward 
shift effect disappears. This observation, which suggests 
a freeze-out curve staying close to the lattice estimate [B[ 
of the parton-hadron coexistence line in the intermediate 
Hb domain, will be further checked as high statistics data 
from the STAR runs during the RHIC beam energy scan 
(BES) program become available. 

An indication of analogous behavior is observed in SM 
analysis of the corresponding SPS data, which needs to be 
confirmed by more detailed analyses, including a study of 
the centrality dependence of the chemical freeze-out tem- 
perature at SPS and the use of the core-corona model to 
fit the data. Also, the observed flatness of the central- 
ity dependence of antibaryon yields per participant, in 
minimum bias collisions at this energy, may point to the 
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presence of such a selective suppression. This effect was 
underestimated in our previous publication [3ljj , leading 
to a premature conclusion concerning the existence of an 
alternate QCD phase. 

In summary these observations indicate that the 
hadronic freeze-out curve shown in the QCD phase dia- 
gram of Fig.[T]requires revision, at least, in the interval of 
SPS energies, at baryochemical potential from about 250 
to 450 MeV, where the "empirical freeze-out curve" from 
statistical model analysis begins to disentangle from the 
parton-hadron coexistence line proposed by recent lattice 
QCD extrapolations to finite baryochemical potential. 
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